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Abstract 


We  study  the  limiting  behaviour  of  large  systems  of  two  types  of  Brownian  particles 
undergoing  bisexual  branching.  Particles  of  each  type  generate  individuals  of  both  types, 
and  the  respective  branching  law  is  asymptotically  critical  for  the  two-dimensional  system, 
while  being  subcritical  for  each  individual  population. 

The  main  result  of  the  paper  is  that  the  limiting  behaviour  of  suitably  scaled  sums  and 
differences  of  the  two  populations  is  given  by  a  pair  of  measure  and  distribution  valued 
processes  which,  together,  determine  the  limit  behaviours  of  the  individual  populations. 

Our  proofs  are  based  on  the  martingale  problem  approach  to  general  state  space  processes. 
The  fact  that  our  limit  involves  both  measure  and  distribution  valued  processes  requires  the 
development  of  some  new  methodologies  of  independent  interest. 
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Introduction 

In  this  paper  we  shall  study  the  limiting  behaviour  of  large  systems  of  two  populations  of 
Brownian  particles  undergoing  bisexual  branching.  As  an  example,  consider  the  following 
system: 

Assume  that,  at  time  t  =  0,  n  >  1  particles  of  two  types  (type  1  and  type  2)  are  located 
in  and  begin  moving  as  independent  Brownian  motions.  Each  particle  of  type  i,  t  =  1, 2, 
dies,  independently  of  the  others,  after  an  exponential  time  with  mean  1/n.  Assume  that, 
at  the  time  of  death,  each  particle  is  replaced,  on  average,  by  one  particle  of  each  type, 
with  overall  probability  and,  with  probability  by  nothing.  The  replacement  particles,  if 
there  are  any,  perform  independent  Brownian  motions  starting  at  the  point  of  death  of  their 
parents,  and  this  story  of  diffusion  and  (critical)  branching  continues  until,  with  probability 
one,  there  are  no  particles  left. 

Let  •)  be  the  measure  valued  processes  describing  the  positions  of  the  particles  at 

time  i.  That  is,  let 

A)  =  Number  of  particles  of  type  t  in  A  at  time  i. 

Our  interest  lies  in  establishing  the  joint  behaviour  of  the  Zj”\t)  in  the  infinite  density, 
n  — ♦  oo,  limit. 

This  problem,  in  the  absence  of  particle  motion,  has  already  been  studied  by  Kurtz  [5]. 
The  addition  of  the  particle  motion,  however,  makes  the  problem  considerably  more  complex, 
and  rather  interesting. 

In  order  to  describe  a  typical  result,  define  the  following  two,  rescaled,  processes. 

XW(()  =  Z.W(I)  +  zW(l)^  ^  z!">(t)  -  zW(i) 

n  n 

We  shall  study  the  limiting  behaviour  of  and  as  n  oo.  This,  clearly,  will  tell 
us  about  the  limiting  behaviour  of  the  Zj^K  What  we  shall  find  is  that  while  has  an 
interesting  nontrivial  limit,  converges  to  zero.  That  is,  in  the  infinite  density  limit  the 
proportions  of  puticles  of  each  type  are  identical. 

This,  naturally,  leads  one  to  try  to  re-rescale  the  difference  so  as  to  obtain  a  fluc¬ 
tuation  result  describing  the  rate  at  which  balance  between  the  two  populations  is  achieved. 
It  turns  out  that  this  is  best  done  via  a  third  process,  defined  by 

t^(n)(<)  =  y(»)(f)  +  /‘„y(’*)(5)rf5. 

Jo 

We  shall  show  that  Wf")  has  a  nice  limit,  as  n  — »  oo,  as  a  (Schwartz)  distribution  v^ued 
process,  and  that  the  convergence  is  joint  with  that  of  X^^K 

The  main  result  of  the  paper.  Theorem  3.3,  gives  the  details  of  this  convergence,  under 
a  more  general  setup  than  that  just  described. 

Since  our  study  involves  both  measure  and  distribution-valued  processes  we  shall,  un¬ 
fortunately,  need  to  start  with  some  technical  results  about  weak  convergence  for  measure 
cross  distribution  valued  processes.  Some  of  these  results  should  be  of  independent  interest. 
This  is  Section  1  of  the  paper. 
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In  Section  2  we  begin  describing  our  system  in  detail,  and  describe  also  the  previous 
result  of  Kurtz  [5]  noted  above.  Section  3  contains  the  main  result  of  the  paper,  and  Section 
4  is  devoted  to  proofs.  There  is  an  appendix  containing  the  proof  of  the  existence  and 
uniqueness  of  the  solution  a  particular  non-linear  evolution  equation  which  is  used  in  the 
proofs  of  Section  4. 

1  Preliminaries,  weak  convergence 

In  this  section  we  give  a  brief  introduction  to  the  meu'tingade  problem  approach  to  general 
state  space  process.  We  start  with  some  notation  and  definitions.  Let  E  be  topological  space, 
B{E)  {C  {E))  be  the  set  of  bounded  (bounded  continuous)  Borel  measurable  functions  on  E. 
Denote  by  B(E)  t!.e  <T-algebra  of  Borel  subsets  of  E  and  by  V  (E)  the  set  of  Borel  probability 
measures  on  E.  Let  A  be  subset  (not  necessary  linear)  of  B{E)  x  B{E). 

Definition  1.1  By  a  solution  of  the  martingale  problem  for  A  we  mean  a  measurable  stochas¬ 
tic  process  X  with  values  in  E  defined  on  some  probability  space  (n,.F,  P),  such  that  for  each 
ill  9)  ^  ^  process 

(1.1)  im))-  f9ms))ds 

Jo 

is  a  martingale  with  respect  to  the  filtration 

a  (^J\  (Xiu))  du:  s<t,  he  B{E)^ 

where  =  <r  (A’(s)  :  s  <t). 

Note  that  if  A"  is  a  right  continuous  process  then  . 

Definition  1.2  When  an  initial  distribution  p  G  "P^E)  is  specified,  we  say  that  a  solu¬ 
tion  of  the  martingale  problem  for  A  is  a  solution  of  the  martingale  problem  for  {A,p)  if 
PA(0)-‘  =  p. 

Let  {.Ft}  be  a  complete  filtration.  Let  C  be  the  space  of  progressive  (i.e.  {Ft}-progressive) 
pnx^sses  Y  such  that  supt>o  E  [|K(t)|]  <  00.  Set 

(1.2)  ||K||  =  sup£;r(0|]. 

f>0 

We  do  not  distinguish  between  C  and  the  Banach  space  of  equivalence  classes  in  £,  deter¬ 
mined  by  the  norm  (1.2).  (  A  ~  K  if  |(A  —  K||  =  0  ).  We  define  the  semigroup  of  operators 
S{s)  on  C  by 

(1.3)  (J(s)r)(0  =  £^[i"(«  +  5)|:Ft]. 

Define 

(1.4)  A  =  |(y,  Z)eCxC:  Y{t)  -  Z{s)  ds  is  Fi-martingale| . 
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A  is  called  the  full  generator  of  Jis). 

Let  ^(£>£[0,  oo))  be  the  a-algebra  generated  by  the  simple  cylindrical  subsets  of  Z?£:[0,  oo); 
i.e. 

g{DE[0,oo))  =  a  ((tt., . {BiE))^  1 e  [0,oo),  m  e  AT)  , 

where  for  points  in  the  projection  ittyM . -DeIOiOo)  — *  £'"*  is  defined 

by 

^tuu . «m(^)  =  W<i),ar(t2). •••,!(<»»)),  Vx  €  Z)e[0,oo). 

When 

(1.5)  a(f?E(0,oo))  =  5(Z?£[0,oo)) 

we  shall  say  that  E  has  equivalent  Borel  and  cylindrical  a-algebras. 

We  now  give  an  analogue  of  the  Theorem  4.8.10  [5]  without  the  assumption  that  the 
space  E  is  metric.  The  result  will  be  crucial  for  our  later  needs,  as  our  main  result  will, 
ultimately,  follow  directly  by  checking  that  the  conditions  of  this  result  are  satisfied. 

Theorem  1.3  Let  (E,t)  be  a  completely  regular  topological  space  with  equivalent  Borel  and 
cylinder  c -algebras  on  D£;[0, oo).  For  each  n  >  \  let  be  a  complete  filtration,  and 

let  A  C  ^ (E)  X  ^ (E)  and  v  €  'P{E).  Suppose  X„,  n  =  1,2, ...,  is  a  -adapted  process 
with  sample  paths  in  oo),  {A’„}  is  relatively  compact,  and  PXn{0)~^  =>  v,  as  n  oo. 
Suppose  furthermore  that  for  each  (/, h)  €  A  and  T  >  0,  there  exist  (^n>‘Pn)  €  A„,  (where 
An  is  defined  as  at  (l.f)  but  for  Tf  martingales),  such  that  the  following  three  conditions 
hold: 

(1.6)  supsup£;[|^„(s)|]  <  00, 

n  s<T 

(1.7)  supsup£;[|<^„(s)|]  <  OO, 

n  s<T 

(1.8)  lirnE  [\Us)  -  f  {X„m]  =  limf;  I|^„(s)  -  h  (X„(0)|]  =  0. 

nioo  njoo 

Then 

(a)  Each  limit  point  of  {Xn}  is  a  solution  of  the  f?£;[0,oo)  martingale  problem  for  (A,!/). 

(b)  If,  in  addition,  we  assume  that  the  2?e(0,  oo)  martingale  problem  for  (A,  u)  has  at  most 
one  solution,  then  Xn  =>  X ,  as  n  —*  oo,  where  X  is  the  unique  solution  in  De^i  oo) 
of  the  martingale  problem  for(A,v). 

Proof  The  proof  that  each  limit  point  of  {A'„}  is  a  solution  of  the  martingale  problem  for 
{A,t/)  will  be  analogous  to  the  proof  of  Theorem  4.8.10  [5]. 

Let  K  be  a  limit  point  of  Let  (/,A)  €  A  and  T  >  0,  and  let  ((n,Vn)  satisfy  condi¬ 

tions  (1.6  -  1.8).  Let  t>0,  0<ti<t2<...<tk<t<t  +  s<T.  Since  (^„,<,s'n)  €  A„,  it 
follows  that 

(1.9)  £[(«.(<+>) -WO -/%.(»)  =»• 
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for  all  ft}, . . . ,  Afc  €  B{E),  n^^f,  which  implies  that 

£  +  «)  -  fiXnit  +  «))  -  (^n(0  -  fiXnit)))  ~  /%»(«)  -  A(Jr„(«))  rfti]  (Jf„(t.)) 

+  E  |(/(Xn(<  +  3))  -  f{Xn{t))  -  J**\iXn{u))  du)  j[hi  (X„(t.))]  =  0,  Vn. 

Taking  into  account  condition  (1-8),  and  that  ft}, . . . ,  ft^  6  B{E),  we  have  that 

(1.10)  Hmf;  J(/  {Xnit  +  s))  -  fiXnii))  -  (Xniu))  du)  j[hi  (XniU))  =  0 
and  hence 

(1.11)  B[(/(n<+»))-/(i'(<))-/*'A(i'W)<i«)n*.(nii))  =0. 

Thus  we  have  that  Y  is  a  solution  of  the  mairtingale  problem  for  (4,  i/),  and  part  (a)  is 
proved. 

By  (a)  and  (b)  we  have  that  the  martingale  problem  for  (^4,  i/)  has  a  unique  solution. 
This  means  that  the  finite-dimensional  distributions  of  all  limit  points  of  {Xn}  coincide. 
Since  the  Borel  and  cylinder  o'-algebras  are  equivalent,  the  finite-dimensional  distributions 
determine  the  probability  measure  on  DglO,  oo).  By  this  argument,  together  with  relative 
compactness,  we  get  that  each  subsequence  of  {A’n}  contains  a  further  subsequence  which 
converges  to  the  unique  solution  of  the  martingale  problem  for  (A,v).  By  Theorem  2.3  [2] 
we  are  done.  ■ 

In  order  to  apply  (b)  of  Theorem  1.3  we  need  to  know  how  to  determine  uniqueness  for 
solutions  of  martingale  problems  for  processes  with  values  in  ubitrary  topological  spaces. 

Theorem  1.4  Let  (E,t)  be  an  arbitrary  topological  space,  and  let  A  C  B(E)  x  B(E). 
Suppose  that  for  each  fi  {E)  any  two  solutions  X,  Y  of  the  martingale  problem  for 
{A,  ft)  have  the  same  one-dimensional  distributions;  i.e.  for  each  t>  0  , 

(1.12)  p{X(0€r}  =  p{Y(t)€r},  r€B(£:). 

7%en  any  two  solutions  of  the  martingale  problem  for  {A,  p)  have  the  same  finite-dimensional 
distributions;  i.e.  (1.12)  implies  uniqueness  on  the  cylinder  a -algebra.  If  X  is  a  solution  of 
the  martingale  problem  for  {A,  p)  with  respect  to  the  filtration  Tt ,  tften 

(1.13)  B  U  (X(,  +  ())  l^-J  =  £  1/  (X(5  +  ())  |.)f  (.)) 
for  all  f  £  B(E)  and  s,t>0. 

Proof  The  proof  is  completely  analogous  to  the  proof  of  Theorem  4.4.2  of  [5]. 

Corollary  1.5  Let  (E,  r)  be  a  completely  regular  topological  space  with  equivalent  Borel  and 
cylinder  a-algebras  on  j!?E[0,oo).  Let  A  £  B{E)  x  B{E).  Suppose  that  for  each  p  6  ViE) 
any  two  solutions  X,  Y  of  the  martingale  problem  for  {A,p)  with  sample  paths  in  Z7£[0,oo) 
satisfy  (1.12)  for  each  t  >  0.  Tften,  for  each  p  £  V{E),  any  two  solutions  of  the  martingale 
problem  for{A,p)  with  sample  paths  in  f?£[0,oo)  have  the  same  distributions  on  P£[0,oo). 
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Proof  The  fact  that  X  and  Y  have  the  same  finite-dimensional  distributions  is  a  part  of 
Theorem  1.4.  Since  the  Borel  and  cylinder  tr-algebras  are  equivalent  the  finite-dimensional 
distributions  of  X  and  Y  determine  their  distributions  on  Z?£;[0,  oo),  and  we  are  done.  g 

Let  denote  the  space  of  hnite  measures  on  (^,  endowed  with  the  topology 

of  weak  convergence.  For  n  €  and  /  €  let 

(114)  (/,(<)  2  jidii. 

Denote  by  S{Rf)  the  (Schwartz)  space  of  rapidly  decreasing  functions  on  i?**,  and  by 
S\R^)  the  topological  dual  of  the  space  of  tempered  distributions.  We  endow  S'{B^) 

with  the  strong  topology. 

For  reasons  that  will  become  clear  later  we  now  need  to  study  stochastic  processes  taking 
values  in  A/f  (/?'*)  x  S\R^).  For  simplicity  denote  Mp(R^)  x  by  Mf  x  S'.  Recall  that 

in  the  theory  of  weak  convergence  of  processes  in  Z)£[0,  oo),  where  E  is  complete  separable 
metric  space,  the  basic  equivalence  (1.5)  between  the  Borel  and  cylindrical  o-algebras  is  of 
crucial  importance.  It  is  not,  however,  clear  that  it  carries  over  to  the  case  of  E  =  Mf  x  5', 
which  is  what  we  shall  require.  The  general  problem  of  determining  which  spaces  E  do  possess 
equivalent  Borel  and  cylindrical  o-algebras  has  been  studied  in  some  detail  by  Jakubowski  [7]. 
Following  his  techniques,  we  shall  establish 

Theorem  1.6  The  space  DAf^x5'[0,oo)  has  equivalent  Borel  and  cylinder  a -algebras. 

Before  the  proof  of  this  Theorem  we  note  the  following  lemma  from  [7]. 

Lemma  1.7  (i)  Suppose  that  the  completely  regular  topological  space  {E,  r)  has  the  fol¬ 

lowing  two  properties: 

(1)  Compact  subsets  of  E  are  metrizable. 

(2)  There  exists  a  sequence  of  of  compact  subsets  of  E  such  that  for  every 

X  €  De[0,  1]  one  can  find  A'„  containing  the  set  x  =  {x(<)|f  €  [0, 1]}. 

Then 

(1-15)  a(X?£l0,l])  =  B(DE[0,l]). 

(ii)  If  a  completely  regular  linear  topological  space  E  satisfies  (1.15)  it  also  satisfies  (1.5). 

Proof  of  Theorem  1.6  First  let  us  prove  that  Mf  x  S'  generates  equivalent  Borel  and 
cylinder  o-algebras  for  xS'[0, 1]. 

We  know  that  Mf  is  a  separable  metric  space,  hence  it  is  homeomorphic  to  a  subset  of  R". 
The  space  R  has  the  properties  (1)  and  (2)  from  Lemma  1.7,  and  by  Proposition  5.3  [7]  S'  also 
satisfies  these  conditions.  Hence  by  Corollary  2.8  [7]  R°°  x  S'  has  the  property  (1.15)  and 
by  Theorem  2.1  [7]  Ax  S'  also  has  this  property  for  each  subset  A  of  By  Jakubowski’s 
Theorem  1.3  the  topology  on  1]  depends  only  on  the  topology  t  on  E,  consequently 
the  property  (1.15)  is  preserved  under  homeomorphism.  Thus,  by  a  homeomorphism  of  Mf 
to  some  subset  of  we  have  that  (1.15)  is  satisfied  by  Mf  x  S'  and  hence  by  part  (ii)  of 
Lemma  1.7  i5jvfj,.x5'[0»«>)  has  equivalent  Borel  and  cylinder  o-algebras,  as  required.  | 

We  need  some  additional  properties  of  the  space  Mf  x  S'. 
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Lemma  1.8  (i) 


B(Mf  X  S')  =  B{Mf)  X  B(S'). 

(ii)  Ltt  N\  C  U(Mf)  and  C  U(S')  be  aeparnting  for  V{Mf)  and  V  (S')  respectively. 
Then 


^  =  {/i/2:/i€MU1,  /,e  ATjUl} 
is  separating  for  V  (Mf  x  S'). 

Proof  (i)  It  is  sufficient  to  show  that 

(1.16)  B(Mf  X  S')  C  B{Mf)  X  5(S'). 

Let  be  a  Borel  subset  of  Mf  x  S'.  By  Proposition  5.3  [7]  we  have  that 

(1.17)  A=  \J  Ar\{MFxKn), 

where  Kn  are  compact  in  S'  and  metrizable.  Each  compact  metric  space  is  separable,  and 
Mf  is  a  separable  metric  space,  hence  by  [2],  p.225  we  have  that 

(1.18)  B  (Mf  X  Kn)  =  B(Mf)  x  B  (Kn) ,  Vn. 

From  the  definition  of  the  relative  topology  in  Mf  x  Kn  we  have 

(1.19)  B(MFxKn)  =  tr(Bn(MFxKn),  Ber) 


=  {Bn(MFxKn),  BeB(MFxS')}, 

where  r  is  the  topology  on  Mf  x  S',  and  the  last  equality  follows  by  [3],  Theorem  10.1.  By 
choice,  A  €  B(Mf  x  S'),  and  so 


Thus,  by  (1.18), 


A  n  (Mf  X  Kn)  €  B  (Mf  x  Kn) . 
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An(MFX  Kn)  €  B(Mf)  X  B  (Kn) . 
By  the  same  arguments  as  in  (1.19)  we  get 

B(Kn)  =  {BnKn.B^B(S')}, 

so  that 

B(Kn)cB(S’),  Vn, 

and 
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An(MFX  Kn)  €  B(Mf)  X  B(S'),  Vn. 


By  (1.17)  we  obtain 


A  e  B(Mf)  X  B(S'), 


and  we  are  done. 

(ii)  The  proof  of  this  part  is  completely  analogous  to  that  of  Proposition  3.4.6  [5].  g 

Let  Ci{R^)  denote  the  set  of  continuous  functions  with  limit  at  infinity.  In  general,  if  F 
is  a  set  of  functions  on  if*,  write  F+  for  {/  e  F  :  infjjd  /(x)  >  0}. 

Corollary  1.9  The  set  of  functions 

{Z'  G  7f(MF  X  S') :  />,./,  =  exp  {(/,,//,)  +  i  (/2,a<2)}}  i 

where  f\  G  and  /j  G  S(R^),  is  separating  on  P(Mf  x  S'). 

Proof  By  [4],  Theorem  3.2.6,  the  set  of  functions  {exp{(/,  •)},  /  G  C/(ff*)+}  is  separating 
on  P  (Mf)-  By  [6],  Theorem  3.2,  the  probability  law  of  X,  where  A’  is  a  random  distribution, 
is  uniquely  determined  by  the  characteristic  functional 

Cx(/)  =  £[exp{t(/,A)}],/G5(/?-'), 

and  hence  the  set  of  functions  exp  {t  (/j,  •)}  is  separating  on  P  (S').  The  result  then  follows 
by  Lemma  1.8,  in  spite  of  the  fact  that  1  ^  {exp {(/,  )},  /  G  since  there  exists 

{/n}  €  C/(/2^)+,  such  that 

(1.20)  bp-^Iijn  exp {(/,  )}  =  1. 


2  Bisexual  branching  without  diffusion 

In  this  section  we  shall,  briefly,  describe  some  results  of  Kurtz  [5],  which  correspond  to  a 
generalization  of  the  bisexual  branching  system  described  in  the  Introduction,  but  for  which 
the  “particles"  perform  no  motion.  In  the  following  section  we  shall  extend  this  model  to 
the  one  of  interest  to  us,  but  it  is  worthwhile,  at  this  stage,  to  look  at  the  simpler  case. 

Kurtz  [5]  considered  a  system  made  up  of  two  types  of  particles.  E^ach  particle  lives  for 
an  exponentially  distributed  lifetime  with  parameter  Aj  or  A2,  depending  on  its  type.  If  a 
t3rpe  1  particle  dies  it  gives  rise  to  offspring  of  types  1  and  2  with  the  number  of  offspring 
distributed  as  (71,72).  Similarly,  if  a  type  2  particle  dies,  it  gives  rise  to  offspring  of  both 
types  distributed  as  (V’i>V'2)-  Assume  that  J5[7j®]  <  00,  £[0f]  <00,  t  =  1,2.  Define 

^  bj] .  rntj  =  E[i>j],  j  =  1 , 2. 

Let  Zj  be  the  number  of  type  j  particles  alive  at  time  t,  and  set 
(2.1)  =  t>B. 
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Then  Z  is  is  a  two-type  Markov  branching  process.  We  shsdl  assume  that  the  process 
is  critical  and  rriij  >  0  for  all  t,j  =  1,2.  This  implies  that  there  exist  vectors  (i/i,!/])  and 
and  a  real  number  ij  >  0  satisfying 


(2.2) 

(  Ai(mii-l) 

Aimi2 

^  A2m2i 

3 

1 

(2.3) 

/  Ai(mn  -  1) 

Aimi2 

V  A2m2i 

^2(^22  ~  1) 

Taking  i/i,i/2  >  0  and  and  ^2  will  have  opposite  signs.  Kurtz  considers  a  sequence  of 
processes  ,  Z^"^)}  with  initial  population  sizes  (Zj  (0),  Z2^\o))  =  ([nzi],  [7122]),  where 

(21,22)  are  fixed.  Set 


(2.4) 

(2.5) 


A'W(t)  = 

n 


Then  both  and  K^’*^(<)  exp{njj<}  are  martingales  [1].  We  are  interested  in  the  limiting 

behaviour  of  and  as  n  — ►  00.  As  might  be  expected  from  the  fact  that 

K^"^(f)exp{nf;t}  is  a  martingale,  ^^"^(f)  converges  to  zero  as  n  — >  00,  and  so 


2!'"("t),.  6XW(<) 

n 


_ ^  zr\nt)  irX^-){t) 

1^1(2  -  ‘'26  ’  n  1/2(1  -  ui(2 


Thus  we  have  the  rather  surprising  result  that  the  limiting  behaviour  of  gives  the 
limiting  behaviour  of  both  Zi'\nt)/n  and  Z^\nt)/n  for  t  >  0. 

The  limiting  behaviour  of  is  somewhat  more  delicate,  and  is  best  described  in  terms 
of  a  new  process  defined  by 

Wt">(0  =  Kl">(t)  +  /‘niyyt">(s)  ds. 

Jo 

To  state  Kurtz’s  main  theorem  we  require  the  random  variables 

7i  =  *'1(71  -  1)  +  *'272  ,  72  =  6(71  -  1)  +  672, 

(2.6) 

^1  =  I/iV*!  +  V2{^2  -  1)  ,  ^2  =  Ml  +  6(02  -  !)• 


and  the  parameters  ajj  =  and  a?-  =  £[^,-0j].  Then  Kurtz  [8]  proved 

Theorem  2.1  (a)  The  sequence  {(X^’*\  W^"^)}  converges  in  distribution  to  a  R^-valued 
diffusion  (X,  W)  with  generator 

X 

(2.7)  Af{x,  tn)  =  -  (an/t*(i.  «*)  +  2aj2/*,«,(x,  w)  -I-  a2if^{x,  w)) 
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whert  Oij  =  (Ai^jaJ^  -  A26a?j)/(‘'i6  - 
(b)  For  each  T  >  0, 

sup|K^”^(<)  —  K^”^(0)exp{— ru7<}l  — ♦  0,  in  probability,  asn  —*  oo. 
t<T 

Consequently,  for  0  <  <  tj,  it  follows  that  //*  nqY^'^^s)  ds  converges  in  distribution  to 

W^(^2)  ~  as  n  —*  oo. 

Our  aim  now  is  to  add  diffusion  into  the  above  model,  and  obtain  a  version  of  Theorem 
2.1  incorporating  the  spatial  motion. 

3  Bisexual  branching  with  diffusion 

We  start  with  a  “typical”  situation.  Assume  that  at  time  f  =  0  we  have  a  number  of  particles 
of  types  1  and  2  scattered  throughout  Bf,  which  begin  diffusing  as  independent  Brownian 
motions  with  generators  jA.  Each  particle  of  type  t,  t  =  1,2,  dies,  independently  of  the 
others,  after  an  exponential  time  with  parameter  A,.  Let  be  the  probability  that  at  death 
a  particle  of  type  i  produces  k  offspring  of  type  1  and  /  offspring  of  type  2.  For  Borel  A  C 
set 


Zi{t,  A)  =  Number  of  particles  of  type  i  in  the  set  A  at  time  t. 

Then  Z,(t)  €  Mp  for  all  <  >  0.  Consider  the  vector  measure-valued  process 

z(t)  =  {^,(0,Z2(0), 

adapted  to  the  filtration 

Tt.  =  a  ^Zs  :  s  <  tj  . 

We  write 

(3.1)  \Zi{t)\  =  i  =  l,2, 

to  denote  the  total  mass  processes,  essentially  equivalent  to  the  population  processes  dis¬ 
cussed  in  the  previous  section. 

Take  now  a  sequence  Z^'^\t)  =  {Zi'\t),  Z^\t))  of  such  processes  with  death  inten¬ 
sities  Ai"\  A^"^,  and  offspring  distributions  pl/"\  adapted  to  the  filtrations  = 

a  :  s  <  t). 

Let  (7i^"^  >72^"^)^  have  joint  distribution  p|y"\  and  (^{"^  have  joint  distribution 

p*j"\  Set 

mW  =  £[7/"l),  ra<”>  =  £(*(■>]. 
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Let  denote  eigenvectors  corresponding  to  the  eigenvalues  , 

—172“^  of  the  matrix  where 


(3.2) 


A<">  = 


(  - 1)  Ai"’mi»2>  \ 

-  1) 


Then  we  make  the  following  regularity  assumptions: 
,(»)  ^  1  n  ^  ^ 


(3.3)  0<m|"^<l,  Vn,  0  <  =  lim  m,-"^  <  1,  t,j  =  1,2, 

^  fif  00 


(3.4)  aJ"^  >  0,  t  =  1,2,  Vn,  either  lim-f-r  or  lim-7-7  exists,  limA-"'  =  00,  t  =  1,2, 

nlooJ^j”*  nfoo 


(3.5)  lim  -  l)  ~  0  “ 


(3.6)  lim 

’  nToo 


Ai"^a1"^  -  1)  (m^2^  -  0  - 


Ai">(l  -  ml;))  +  Al»)(l  -  ml;)) 


<  00, 


(3.7)  sup  £[(7/"))^  <00,  sup£ 


<  00, 


a(”) 

(3.8)  sup  — —  <00,  i  =  1,2. 

n  n 


Simple  but  tedious  calculations  show  that 


(3.9) 

II 

1 

3 

(3.10) 

’ll"’ 

=  A<”'(l-mW) 

Set 

(3.11) 

f/i  =  limj?!"). 
nfoo 


Under  conditions  (3.3)  -  (3.6)  we  obtdn  that  |f7i|  <  00. 

Despite  the  heavy  notation  above,  things  are  are  not  as  difficult  as  they  seem.  For 
example,  if  we  return  to  the  example  of  the  Introduction,  then  we  have  Aj")  =  Aj")  =  n 
and  m|“)  =  1/2,  *,i  =  1,2.  It  then  follows  that  17I")  =  0,  i/l")  =  n  and  we  can  t«dce 

i/i")  =  1/2")  =  1,  and  ^1")  =  —1,  ^2"^  =  1  for  eadi  h  >  0. 

The  following  technical  lemma  follows  via  straightforward  algebra. 
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Lemma  3.1  Under  conditions  (3.3)  -  (3.6)  there  exists  a  finite  N  >  1  and  eigenvectors 

such  that 

(3.12)  >“{,*'.•"^>0. 

n>N 

(3.13)  <  0.  >0,  V  n  >  JV,  lim  <  oo, 

ntoo  I  I 

(3.14)  lim  -  4"’{ii  >  0. 

Remark  3.2  Without  loss  of  generality  we  shall,  henceforth,  assume  that  =  1  in  (3.12). 


As  in  the  case  of  Kurtz’s  simpler  model,  we  need  to  define  a  few  more  additional  random 
variables  and  coefficients  before  we  can  state  our  main  result: 


(3.15) 


7"  =  -  1)  +  .  7i  s  -  1)  +  ^”^2*”’, 

«  =  +  .'S"’(4"’  - 1) .  «  s  f +  {'"|(^S"|  - 1). 


(3.16) 


(.1  _  -  AWfl-'afl 

"(.i-'d-’-.'rd"') 


a,j  _ 


[  af^  =  E  ,  af^  =  E  ,  »,J  =  1,2. 


Assume  that  limoj”^  <  oo  exists  (this  is  not  a  strong  restriction  given  the  results  of  the 


ntoo 


previous  lemma  and  assumptions  (3.3)  -  (  3.8))  and  define 


(3.17) 


an  =  lima.j 

nfoo  ^ 


(") 


We  are  finally  in  a  position  to  define  the  three  measure  valued  processes  that  interest  us. 
(Note  that  the  second  two  of  these  are  signed  measures,  or  distributions.) 

n 


(3.18) 


'  n  ’ 

WW(f)  =  K(”)(f)  4- jf‘»7fV(")(s)ds  , 


In  general,  let  V{A)  denote  the  domain  of  an  operator  A.  Here  is  the  main  result  of  this 
paper: 
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Theorem  3.3  Let  {(X^"^(0),  W^"^(0))}  have  a  limiting  distribution  v  €  T\Mf  x  5),  and 
assume  that  {|A^^"^(0)|*}  is  uniformly  integrable.  Suppose  that  (S.S)  -  (  S.8)  hold  and 
that  liinn_oo  “Ij  ^  exists  a.id  is  finite.  Then  =>  (A^,  W'),  where  (A,  e 

Cj»/rx5'[05®o)  ^  unique  solution  in  i)A#irxS'[0»oo)  of  the  following  martingale  problem 
for  (A,  v): 

A=  {exp{-(5i,/ii) +  t{52,/^2)}, 

(3.19)  exp{-  {gi,pi)  +  i  {S2,A‘2>}  (5  (-Afli  -  2r)igi  +  a^g^  -  2iaiigxg2  -  a22fl'|)  ,/^i)  : 
SieP(iA)+,fl2  6S(/2‘')}, 

where  Oij  and  t)i  are  as  in  (3.17)  and  (3.11)  and  V{^A)+  =  T>(^A)  fl  Ci{R^)+  . 

Furthermore,  for  each  T  >  0  and  g^  €  S{Ef), 

(3.20)  sup|^52,5^^"n0)  -  (52,V'(’*)(0))exp|-7/^”^i}|  -+  0,  tn  probability, 
as  n  00. 

The  rest  of  this  paper  is  devoted  to  the  proof  of  above  theorem.  The  proof  will  rely  on 
checking  the  conditions  of  Theorem  1.3. 

4  Proofs 

There  are  a  number  of  steps  involved  in  applying  Theorem  1.3  in  order  to  prove  Theo¬ 
rem  3.3.  The  first  of  these  lies  in  finding  the  approximations  and  ipn  to  /(A^"\iy^"^) 
and  ^(At"),  ly^")),  where  {f,h)  €  A.  We  do  this  in  the  following  subsection,  denoting  the 
approximations,  in  Lemma  4.1,  by  and 

Section  4.2  contains  a  sequence  of  preparatory  lemmas  that  ultimately  show  that  (A^"^  W<”^) 
satisfies  a  compact  containment  condition,  required  in  Section  4.3,  for  showing  that  ( A^"\  W^”)) 
is  relatively  compact  in  f?MFxS'[0>oo),  a  necessary  requirement  of  Theorem  1.3. 

The  final  requirement  of  Theorem  1.3,  that  the  martingale  problem  (/I,  v)  have  a  unique 
solution,  is  established  in  Section  4.5  .  This  proof,  itself,  relies  on  a  uniqueness  result  for  a 
particular  non-linear  evolution  equation,  whose  proof  is  relegated  to  ein  Appendix. 

In  Section  4.4  we  show  that  the  limit  process,  which,  according  to  Theorem  3.3  and  the 
previous  steps  is  in  Z)Afpx5'[0>oo)?  is  actually  continuous,  and  that  (3.20)  holds. 

4.1  The  martingale  approximation 

We  need  some  simple  notation  which  will  be  used  later.  If  for  sequence  of  random  variables 
{fn}  there  exists  a  constant  C  >  0,  independent  of  n,  and  another  sequence  {t^n}  such  that 

(4.1)  Knl<C|V>„l,  Vn 

we  shall  write  that  =  O  («/’n)-  Denote  by  o(l)  a  sequence  of  uniformly  bounded  random 
variables  which  converges  to  zero  as  n  — »  00. 


Define  Z  by: 

—  f  f(x  10^  1 

(4.2)  L  =  = 

For  simplicity  we  shall  denote  by  /^ ,  /«, ,  /*, ,  ,  Atu  the  first  and  second  order  partial 

derivatives  of  /  e  Z. 

For  each  /  €  Z,  pi  €  (  where  ^{R^)  is  the  set  of  bounded  continuous  functions 

on  R^  ),  and  gt  €  S{R^)  define 

(4.3)  :MfxS'~^  R 
by 

(4.4)  Gf,g,,g,{fii,  P2)  =  /((fli,Pi),  {92,1^2)),  V(pi,  P2)  e  Mr  X  5'. 

Furthermore,  define  the  operator  Ai  €  Z?(Mf  x  5')  x  V{Mf  x  5')  with  domciin 

(4.5)  I){/1,)  =  :  f€l,g,ev  (ia)^  ,  j,  e  s{k‘)] 

by 

(4.6)  AiGj,g,,gj(fii,fl2)  =  /r((Pl,Pi),(^2,^i2))(»7l<7l,/fl) 

4"  jZ*®  ((5i>  ^1) ,  {52,  A'2))  (an  (Pi)* ,  Pi) 

4“  Zrw  ((PnPi) .  (P2,P2))  (ai2PiP2,Pi) 

4"  jZvnif  ({Pn  Pi> ,  (P2,  P2))  (022  (P2)*  ,  Pi) 

4-  /r  ((PnPi) ,  (P2,P2))  (jApl.pi)  , 

where  the  are  defined  by  (3.17). 

Lemma  4.1  For  all  Gf,gi,g^  €  ^(Ai)  there  exists  (/„,  A„)  €  An,  suck  that : 

m  =  e,,„.„(xi’>,iv<’>}+o((,i">y'ixi’>(t)i), 

K(t)  =  .4,4?,,„,„(A-<">,«'l">)  +  0(n-‘|A'l">(()l)  +  o((i)5">)‘'|.?'">(i)|=) 

+  0((,5">)‘'|X<"'(()|)+0(1). 
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Proof  Set 


/'">(<)  =  /  ((si.  (ft.  »*'<'■>(<)))  ■ 

To  find  80  that  €  -4* ,  we  calculate  limt_oe“‘^ 

obtain 

(4.7)  i<"l(t)  =  (iF'W  ((s.,X("l(l)),(j„W'<")(())) 

+  A<"'£^.  [/  ((s„  XW(0)  +  ^ft,  (ft,  W^'Ht))  +  ^jj) 

-  /((ft.XW(<)).(ft.M'W(()»  ],zi">(o) 

+  {  5^”'"'  {(ft .  X<"'(<)>,  (ft,  H"">(1))) 

+  [/  f(,„XW(())  +  ^ft,  (ft,  +  ^ft) 

-  /  ((ft,  XI“>(1)),  (ft,  ] ,  4'\t) ) 

+  U  ((s., XW(0),  (ft,  H'I’HO))  (7S"’ft, l''"’(0)  "-s- 

where 

(4.8)  (i,  uj)  (•)  =  /r*(x,  tn)  (-|^j 

+  2/xw  (x, »")  ■  ^ digi{-)dig2{  ) 

(n) 

+  /.(x,  u;)^A^,(  ) 

n 

f(n) 

+  /«(x,  u;)  ^Ap,(  ), 

(4.9)  (x,  «j)  (•)  /„(x, 

+  2/,u»  (x,  ti;)  ^ 3,^i(  )9.5j(  ) 

”  <=i 
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+  fww  {x,  w)  S  (%2(  ))^ 

+  /»(*,  ) 

«(") 

+  /„(i,  tr)~Aff2(  ). 

n 

Expanding  /  in  a  Taylor  series  about  have 

(4.10)  A<»>(0  = 

+  ^  [a  ((ji. -«■'"’(<)).  (jj.  »''”’(<)))  E  Hri  SI 

+  A  ((si,xW(i)).(s.,»''">(0))  BRJls! 

+  iA,((si.X'"'(f)>,(si,H''"'(f)»  E  [(7r)1  (S,f 
+  ((si.X‘">(()>.(sj.  *»''■*(<)))  £  [(%")’]  (Si)’ 

+  i/..  ((si.  xW(0),  (si,  »''">(<)))  WiJlsiSi] ,  2l"’(0) 

+  {kE»->  {{gi,X(">(>Ms2,»'^’'it))) 

+  ^  [f.  ((si,X'">(<)).(s!.W'‘">(l)»  E  («J  S. 

+  U  ((si,  A-<"l(()),(s!,  «"">(<)))  «  [«]  Si 
+  ((s..,>f<">(<)>.(si,v»''"'(0))  E  [(«)’]  (s.)’ 

+  ^A.  ((si.  A-">(<)>,  (si.  »''">(<)))  E  [(«)’]  is,)’ 

+  iA.  ((si.X<“l(<)>,(si,H'<”>(<)»  C  [««]  SiSi] .  4”’(o) 

+  A  ((si,  Jf ‘"'(0),  (si.  (sS"'si,l'‘”(<))  +  n-'O  (|X<”'(I)|)  . 
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The  error  is  0(n“*|A’^'‘^(t)|),  since  n~*|Z,-"^(<)|  is  bounded  by  a  constant  times 
Recalling  that 


(4.11) 

(4.12) 


4")(t) 


n 


n 


1 


y 


(4.13) 


we  see  by  (4.13)  that 


(n)  (n) 
vi  m 

Ai"' 


^■K)  =  - 


Ai"> 


~w 


(4.14)  m  52,  Z;"^(t)^  +  [0?]  52,  Y^-\t)) . 

Hence  the  terms  involving  fu,  cancel  and  (4.10)  can  be  rewritten  as 
U^\t)  = 

1  /  "  («S*’f  ■<”)  ((j„  Jf (”1(0),  {»,  iv'">(i)))  -  {|"'F'<">  ((j. , Xl"'(( )),  (j.,  H'W(()» ) 

2'  i,|”>£(">  -  i/("'(j“> 

+  2/.  ((ji ,  -<'("'(1)),  {32,  «''"’(<)» 


+  f..  «9i,X(")(0),(9„H''"'(0» 

(4.15)  +  2/.„  ((j, ,  Jf (“'(()),  {j„  H'(")(0))  oi;>j, j, 

+  U  ({». ,  A^'">(0>.  (S2,  H'("1(0))  «&’  (!»)’ .  X'’>(0  > 


+  /«  ((9i,2f'"'(0).(92.'»^‘"’(<)»  («i)' 

+  2/„  ((ji,  Jf'"'(0),{9..  "^"(O))  ‘IVsift 

+  /w.  ({9.,  X(")(0>,  (9,.  »''">(<)))  tS’  (9.)’ ,  )  +  O  (n-'|2f'"l(t)|) , 
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where 


(4.16) 


and  oj"\  are  defined  in  (3.16).  Note  that,  by  3.8  and  Lenuna  3.1,  the  6|"^  are  uniformly 
bounded  in  n.  Finally,  by  (4.8)  and  (4.9),  we  obtun 

(4.17)*<“>(i)  =  i (/.  ({«„  A-WCi)), (ft,  irl">(l)))  Aft 


+  2A  ((ft ,  XC)(()),  (ft,  w'-ico))  >i!*’ft 


+  ((ft,  («!,  h'I’HO))  <■!?  (ft)’ 

+  2/™  ((ft,2:("l(()>,(ft,H'l">(()))  «!;>ftft 
+  /w  ((ft ,  X(”)(l)),  (ft,  1»'<"'(1)))  (ft)’  ,  Jf'"'(f  )  ) 

+  j  (  /»  ((ft,  X<")(0),  (ft,  WM{t)))  Aft 
+ ((ft,  If  <“>(<)),  (ft,  «''”>(<)»  ‘ff  (ft)’ 


+  2/.„  ((ft,  X<")(l)),  (ft,  H'W(!)))  6';>ftft 

+  U  ((ft,  2f<“>(l)).  (ft,  «''">(!)))  (ft)’ ,  >''"'(<) )  +  O  (»-■  |X(")(()|)  . 

By  (4.17)  it  is  easy  to  see  that 

(4.18)  =  O  (|X<">(0|)  . 

Set 

J<’‘>(x,  U))(-)  =  I  (/«  (x,  to)  Aff2(-)  +  /„  (x,  to)  (ffi(  ))* 

+  2/«w  (x,  to)  ^2^ffl(')j2(*)  +  /fmi  (®»  4j^52(  ))  , 

and 

(4.19)  =  + 

To  find  A<")(<)  so  that  (/^"^  A<"1)  €  A.  we  calculate  liint_K)C"*£^I/^"H^  +  «)  - 
and  obUun 

(4.20)  *<">(0  =  h<")(/)  +  0  !  Y<")(t)l’)  +  O 

-  ((ih,XW(f)),  (52,VV(")(0))  ,  ^^"’(t)) . 
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Combining  tenns  we  obtain  the  desired  result. 

Conditions  (1.6)  (1.8)  now  follow  immediately  from  the  following  bound. 

Lemma  4.2  £[sup,<7.  ts  uniformly  bounded  in  n  for  each  T  >  0. 

Proof  FVom  (4.7)  we  have  that,  for  all  e  >  0, 

(4.21)  Mi*)(0  =  exp{-£|A’(")(0|*} 

+  Jexp  |-£  (5)|  +  I  -  exp  {-£  (s)|*}j  ,  4")(5)^  ds 

is  a  martingale.  Hence 

£[£|.yW(0f  exp  {-£|x(’‘)(0f}] 

<  £[l-exp{-£|j:(")(0f}] 

=  l-£[exp{-£|x(")(0)f}] 

-  E  [j['  [exp  |-.  (|ArM  (,)|  +  |:)’|  -  exp  {-<  wf  }| , 

+  [“P  (IJT*"'  (»)|  +  ^ j  I  -  exp  |jrl">  (.)|’| j  ,  zi%)j  *1 

=  I-E[exp{-.|x'"l(0)|’}] 

+  £[/;(Ai«'  exp  {-£|.^f(")(s)|’} 

X  Jl  -  exp  |-e  (s)|  +  ^)*  -  («)f)  |]  ,  z{"’(«)^ 

+  ^A^">exp|-€|jr(")  (s)|’J 

l-exp|-£^^|x<-)(s)|  +  ^^  ,Z<">(s)^dsj 

<  l-.E[exp|-f|A^<»)(0)|’JJ 
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,(4->«p(-<lx<-'(.)r).^^  [(l*‘*>  wi4)’  -  ^ 

<  c£[lxt->(0)f] 

+.£(/;«p{-<lx“'(.)r}  lW)’ll^l 
+  2^£*.  (*fl  |x<“>(>)|14"’l*)l)  *]■ 

B, ...  ..'.L  «' 

Bdns  a  coMtant  C, ,  indepeadea*  of  a,  aachlk 

(4.22)  [‘  ‘‘'f  {'•  I*'"’  ‘‘'I  J) 

<  .f:[lxi“>(0)f]  +  <c.//ll*‘”’(*)ll''’ 

4  2»,S*> />  [«P  {"  l^'"’  Wf}  ('fl  ‘^• 

H«,maB  .ka.  ,X.-H.)l«P<-W>.  • 

(4.23)  ^  ^  )i 

<  .£  (|XW  (0)f]  4  .<?./;  «.P  R>-}£  O'll  ■" 

+2«,!->j[‘c[«p{-‘i^‘”<'>r}K’‘’>n* 

<  «£  [jxW  (0)f]  +  <C,.«P  ll*'-’  (»)11  . 

By  (4.23)  and  Gronwall’s  inequality  we  obtain 

(4.24)  E[.lxl“)(()f«p{-'1^‘’’(')|  )] 

<  .(£[lX<»(0)f]  4C..«p{.,!’'*}£ll>^‘”(»)ll)“' W"’‘> 


(4.2J) 
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Dividing  by  e  and  letting  c  j  0  gives 

(4.25)  £[|.^(»)(t)|*]  <  (£[|jf(-)(0)f]  +C,texp{i7{")t}£[|A-<»)(0)|])exp{2i7l"’t}. 

That  is,  £[|^(»)(t)|>l  is  uniformly  bounded  in  n .  Once  again  using  the  martingale  properties  of 
we  obtain,  by  Doob’s  inequality,  that  £(sup,<j- is  also  uniformly  bounded  in 
n  for  e^  T  >  0.  m 


4.2  Compact  containment 

Definition  4.3  Let  (E,r)  be  a  completely  regular  topological  space  toith  metrizable  compacts.  Sup¬ 
pose  that  A't"),  n  =  1,2,...,  ts  a  Ep -adapted  process  with  sample  paths  in  D£[0,oo).  We  say  that 
the  compact  containment  condition  holds  for  {Jfl"!},  if,  for  every  e  >  0  and  T  >  0,  there  exists  a 
compact  set  r<,7>  C  £  for  which 

infP{A'(")(0€r..r  for  all  0  <  <  <  t}  >  1  -  f. 

In  the  following  lemmas  we  shall  denote  (dig)^  by  {dg)^. 

Lemma  4.4  ^  bounded  uniformly  in  n  for  each  g2  €  S(Jt^). 

Proof  It  is  easy  to  see  that  for  each  <  >  0 

exp|-« 


FVom  (4.7)  we  obtun  that 

Afp)(t)  =  exp  {-€  (|A'(»)(t)|*  +  (ff2,ir(")(0)’)} 

(4.26) 

-  £  exp  {-c  (s)|’  +  (j2,  } 

+  Ai">£.v-  [«P  +  v)*  +  f  } 

-  exp  {-e  (,)[’  +  (p2,  W^W(t))’)  }]  ,  Z[%)^ 

+  ^cexp  |-c  (|jfW(s)f  +  (ir2,H^(»)(s))’)} 
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+  |exp  |-f  (3)|  +  +  ^(ff2,VV(")(s))  +  j  I 

-  exp  (|jf W  (a)|*  +  (52,  W'("J(0)*)  }]  ,  )) 

-  2€(sa,  W(»)(3))  exp  |-c  (a)|*  +  (pa,  W^(»)(a))*)  }»4’*\»2,y ds 

is  a  martingale.  By  simple  calculations  (as  in  the  previous  lemma)  we  obtain  that 
E  [«(52,W^<")(t))%xp  {-€  +  (ffa.H^(»)(t))*)}] 

<  cf;  [o  (|^(»)  (0)f )]  +  .J*E  [o  (|x(»)  (s)f )]  ds  +  c J^E  [0  (s)|)J  d. 

+  ljf  E  [eexp  |-f  (|a-<")(5)|*  +  (ffa.H^f’‘)(a))*)}(ff2,  ^^<">(5))*]  ds 

=  eBn{t)  +  \f^E  [cexp  {-c  (|x(«)(s)|*  +  (pa,  W'W(s)>')}(p2,lV(»)(s))"]  ds, 

where  the  >  0  are  bounded  uniformly  in  n.  (This  follows  from  the  uniform  boundedness 
of  £lsup,<,  |Ar(")(s)p)).  Hence  there  exists  a  function  B  such  that  B{t)  >  B„(t)  for  each  n.  By 
Gronwall’s  inequality  we  obtain 

E  [(p2,W^<’‘^(t))%xp|-f  +  (p2,1^<'*H0}*)}]  <  B„(0exp{it},  Ve  >  0. 

Letting  c  (  0  we  find,  by  monotone  convergence, 

(4-27)  E  [(p2,W^<'‘H0)*]  <  B„(t)exp{lt}  <  B(t)exp{it}, 


so  that  B[(pa>W^^"H0)*]  is  uniformly  bounded  in  n. 


Lemma  4.5 

(4-28)  (p2,  W^W(O)  -  I  £  (Ap2,  y  i">(5 ))  ds 

is  a  martingale  for  each  n  and  pa  €  S{R^). 

Proof  Note  that  exp{— c|/xi|}sin(e(pa,^a))  €  Vi^Ax).  Define  a  martingale  Afc*\t)  as  in  the 
previous  cases  : 

M^\t)  5  exp  (-C  lAi")  (<)|)  sin  (€(pa,»y<’*>(0)) 

(4.29) 

“  jo  4"^^ 
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+  ^eexp  |-f  |jf<">W|}(  cos 

-  f  sin  (e(fl2,»V(")(s))) 

+  [exp  |-«  (s)|  +  ^)  }  sin  (f(fl2,  W"("Hs))  +  ^52) 

-  exp  {-«  (s)|}  sin  (e(p2,  W^^^Hs)))]  ,  ^{"’(s)) 

+  l^exp  |-«  ^|xt">(s)|  +  I  sin  ^e(j2,  VV<"ns))  +  ^52^ 

-  exp  {-e  |jf ("J  (5)|}  sin  ({(ff2,  W^<"Hs)))] ,  ^3’'^^)) 

+  cexp  |-<  (s)|}cos  (c(p2,W^^"Hs)))»?2"\52’^*"ns))‘^s 
Let 

(4.30)  ~  lim 

=  (ff2,W(-)(t))  -  jf‘  ( [^]  g3,  Zl%)j 
+  [^]  ff2, +  I  (As2,  K<"Hs)) 

=  (S2,  W^(">(t)>  -  hj*  (Aff2,y<"Hs))  ds, 

where  the  last  equality  follows  by  definition  of  The  above  limit  exists  almost  surely.  On  the 
other  hand,  it  is  easy  to  see  from  (4.29)  that 

<  |(52,W^<")(0)|  + ds 

where  Cj  >  0  and  C2  >  0  do  not  depend  on  e.  iVom  the  uniform  boundedness  of  £[|Ar(")(t)|^]  and 
^[(^2*  obtain  the  uniform  integrabiUty  of  |Aft®^(t)/cl,  which  means  that  convergence 

is  in  Lt.  Consequently  is  a  martingale.  ■ 

Lemma  4.8  ^l8up,<j>  (p2,  “  bounded  uniformly  in  n  for  each  §3  €  S{R^). 
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Proof  From  the  previous  lemma  we  have  that 


E  sup(ff2,lVW(t))*  =  E  ^sup  +  ljf‘  (Ap2,y(»)(s)) 

<  2E  [.up  (J/W(0)“]  +  2E  [jupl  (jf'  <As.,yW(.)) 

<  8E  [(Mt=»>(T))*]  +  1 II Alkali"  T'^CE  [sup  (ofj 

<  SE  2(p2,W'<’‘>(T))%i^jr^(Aff2,y(’‘)(s))dsj 

+  illAs2l|*T*C£;|sup|x(’‘)(0fj, 


where  C  is  the  constant  such  that 


|(A<,2,yW(t))f  <  ClIAsrjIl*  |x(»)(0f ,  Vn. 


Finally  we  obtain  that 


E  sup(ff2,py(")(t)}j  <  16f;[(52,W^<"H7’))^]  +  l|As2||*r2C’£  sup|xW(0|"  . 
FVom  Lemmas  4.2  and  4.4  we  obtain  uniform  in  n  boundedness  of  £;(supj<2- 


Lemoia  4.7  For  each  gi  €  P(jA)+,  €  S(R^),  the  processes  {gi,X^”^{t)),  (52. 

i{guX^’^\t)),  (p2>  satisfy  the  compact  containment  condition. 

Proof  The  proof  is  immediate  by  the  uniform  boundedness  of 

E  sup(ff2,W'^’*H0y  .  E  sup(5i,XW(0)* 

.<<7'  '  L<<7'  ' 

and  Chebyshev’s  inequality  in  the  case  of  the  processes  (si,X(")(t)),  (52 >^'’^'*^(0)-  I®  of 

the  process  (52>f^^'*^(0))  result  follows  from  the  fact  that  the  product  of  two 

compact  sets  is  compact  in  the  product  space.  I 
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4.3  Relative  compactness 

We  now  check  the  relative  compactness  of  En  passant,  we  shall  show  that  A’*") 

converges  weakly  to  super- Brownian  motion;  i.e.  the  process  with  values  in  Mf  which  is  the  unique 
solution  of  the  martingale  problem  for  A*  with 

A*  =  |exp{-(si,/i>},  exp{-(fli,^>}(-ihffi  +  jOi,  (fli)’ -  :  y,  €  . 

Recall  that  P(|A)+  =  V(\A)  n  Ci(R<‘)+. 

Lemma  4.8  The  following  sequences  of  processes  are  tight  for  each  gi  €  D(5A)+,  <72  €  S{R'^): 

1)  in  Djp[0,oo), 

2)  {(fli,  A(")(t))}  ,  {(ff2,  ^^(t))}  in  00), 

3)  {(si, A(")(t))  +  (52,  W(-)(t))}  in  Pfl(0, 00). 

Proof  2  and  3  are  simple  corollaries  of  1,  which  we  now  prove.  By  Lemma  4.1  we  obtain  that, 
for  all  G/.ai.jj  €  V{Ai),  we  can  choose  the  processes  and  (/^"\  €  An  such  that 

limE  sup|/<"H0-/((<li.^<"K0).  =  limf:  supO 

=  0. 

The  last  line  follows  from  the  uniform  boundedness  of  fJ[sup,<7’  |A("^(t)|].  Furthermore, 
supf;  sup|/i<’*)(0|  <  supf:  sup|i4i/(^Pi,AW(t)^,  (p2,W(")(t)))| 

-l-supJ?  supO(n“* 
n  [j<T  ' 

+  sup£  s«PO^(»4"^) 

+  supE  sup0((i4’‘>)"'|A(")(t)|2)j  +  sup  0(1) 

<  00, 

where  the  final  line  follows  from  the  uniform  boundedness  of  £[sup«7<  |A(")(t)|^],  and  the  bound¬ 
edness  of  functions  in  the  range  of  the  operator  Ai .  Since  this  holds  for  each  choice  of  gi,  g2,  and 
by  Lemma  4.7  we  have  compact  containment,  applying  Theorems  3.9.1, 3.9.4  [5]  now  completes  the 
proof.  ■ 

Lemma  4.8  and  Mitoma’s  theorem  [10]  immediately  yield 

Lemma  4.0  is  a  tight  sequence  of  processes  m  i?5>[0,oo). 

Lemma  4.10  converges  weakly  to  super-Brownian  motion. 


Proof  We  can  take  ^2  =  0  and  thus  obtain  A  C  U(Mf)  x  where  A  C  Ai  is  defined  as 

follows: 

D(A)  =  linear  span  of  |exp{-  {gi,n)},  g\  eV  (5^)^! 
Aexp{-(ji,/x)}  =  exp{-(5i,/i)}^-T7ifli  + 

By  Theorem  2.4  [14]  and  results  of  Perkins  [12]  the  closure  of  A  generates  a  strongly  continuous 
contraction  semigroup  on  C{Mf)  (with  corresponding  super-Brownian  motion  X)  and  D(A)  is  a 
core  for  A.  By  Theorem  3.2.6  [4]  the  set  of  functions  {exp  {-  strongly  separates  points  in 

Mf-  Thus,  all  the  assumptions  of  Theorem  4.8.2  [5]  are  satisfied  and  we  are  done.  g 

Lemma  4.11  is  a  relatively  compact  sequence  of  processes. 

Proof  We  check  the  conditions  of  Theorem  4.6  [7]  (Jakubowski’s  tightness  criterion). 
and  are  tight.  This  implies  that  the  compact  containment  condition  holds  for  each  of  these 

processes  and  by  the  same  arguments  as  in  Lemma  4.7  it  holds  for  the  pair  Define 

the  family  of  functions  F  :  Mf  x  S'  — ►  R  by 

'fVj.aj(j“i.P2)  =  {51, /ii)  +  (52,M2)  ,  Pi  €  2>(5A)+,  52  €  S{R^). 

This  family  separates  points  in  Mf  x  S'  and  is  closed  under  addition.  By  Lemma  4.8  the  sequence 
{/(X("\  is  tight  for  all  /  €  F.  Jakubowski’s  conditions  are  satisfied  and  the  proof  is 

complete.  ■ 

Theorem  1.3  and  Lemmas  4.11,  4.1  now  imply 

Corollary  4.12  Each  limit  point  of  is  a  solution  of  the  xS'[0>°o)  martingale 

problem  for  (Ai,i/). 

4.4  Continuity  of  the  limit  process 

Throughout  this  subsection  we  shall  denote  by  W  one  of  the  limit  points  of  the  sequence 
and  all  results  are  obtained  for  each  limit  point  of  {W^")).  For  simplicity  of  notation  we  shall 
denote  {g,W{t))  and  {g,X{t))  by  W< (y)  and  Xt{g)  respectively,  while  the  increasing  process  of 
each  martingale  Mt  will  be  {M)f . 

We  shall  treat  only  the  continuity  of  the  process  W) ,  since  the  continuity  of  X  follows  from  the 
fact  that  X  is  super  Brownian  motion.  Using  this  result  we  shall  prove  the  last  part  of  Theorem  3.3, 
viz.  that  ~  converges  in  probability  to  zero  for  each  g%  €  S{R'^). 

Lemma  4.13  For  each  g-i  €  S^R!^),  Wt  {g-i)  is  a  square  integrable  martingale  with  increasing 
process 

(4.31)  {W  (52)),  =  022X,  (p|)  ds. 
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Proof  The  uniform  boundedness  in  n  of  £((W'/"^(S2))*]  (obtained  during  the  proof  of  the  compact 
contsdnment  condition)  immediately  gives  the  existence  of  second  moments  of  Wt  (92)-  We  now  use 
the  same  “trick”  as  in  Lemma  4.5. 

(4.32)  =  exp(-€|Xtl)sin(cH^t(p2)) 

-  exp  {-€  \X,\}  (sin  {iW,  (92))^  \X,\ 

-  cos  ((W,  {g2))ai2X,  (pj)  -  sin  (cW,  (g2))^X,  (pj)) 
-€exp{-€|Jf,|}sin(tW,(p2))»?i  lA",!  ds 

is  a  martingale.  Since 

A/(*)rA 

(4.33)  lim  ‘  ^  =  Wt  (92)  as-, 

the  dominated  convergence  theorem  gives  that  Wt  (92)  is  a  martingale.  Now  note  that 

(4.34)  Mp^(i)  s  exp{-€|A’t|)sin(€W<(52)*} 

-  «exp{-£|X,|}  (csin  (cW,(p2)*)^l^»| 

-  2f  cos  (52)*)  W,  (92)012X1  (92)  +  cos  [tW,  (P2)*) a22X,  [gfj 

-  2iW,  (g^f  sin  (cM^,  (s2)*)a22A^,  (p*) 

-  sin  (cW,  (52)*)»h  ds 

is  a  martingale.  Under  the  same  arguments  as  in  (4.33) 

(2) 

(4.35)  W,  (52)’  -  a22X,  ds  =  ,  a.s. 

is  a  martingale  and  so  (JT  (ff2))(  =  Jo  022X,  (gl)  ds. 

■ 

Lemma  4.14  Assume  that  Wo  €  S'.  Then,  for  each  92  €  S(/Z'*),  Wt  (92)  is  a  continuous  martin¬ 
gale. 

Proof  Define 

(4.36)  Wt^  (92)  =  Wt  (92)  -  A”,  for  all  /r  G  R. 
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Once  again  we  use  the  same  “trick”. 

Afp)(t)  s  exp  {-e  \Xt\ -  {wl<  (p,)’  -  (p,)*)} 

-  jf‘cexp  {-t  \X,\  -  (g^f  -  W^’  (flj)*)  }  \X,\  +  2eW?/‘'  (pa)>  0,2^.  (52) 

+  (2VVf  (ff2)'  -  1)  naaX.  (p»)  -  cih  \X,\)  ds 
is  a  martingale.  Under  the  same  cjguments  as  in  (4.33)  we  obtain  that 
A/(3)(t)  =  exp{-(wl<(g2)^-Wi^igi)^)} 

-  jT'exp  {-  {W^  (52)'  -  (92?)}  (2W,^  (92?  -  l)  a22X.  (5I)  ds 

=  liinA/i^^(<),  a.s. 

is  a  martingale.  By  Corollary  2.3.3  from  [5]  we  have  that  the  process 

=  exp  (92?  -  (92?)  -  £  {2W!?  (92?  -  1)  a22X.  (gl)  dsj 

=  exp  {-  (Wf^'  (92?  -  Wo^  (g2?)  +  j^a22X,  (9I)  ds  -  2jf‘a22H^f  (92?  X,  (pj')  ds^ 
is,  at  least,  local  martingale.  By  setting  92  =  m^92  in  the  above,  we  see  that,  for  every  m  £  R, 
(4.37)  exp  |-m*  (52)'  -  (52)')  -  ^22X1  (pf)  dsj 

is  a  local  martingale.  By  Lemma  4.13 

W!^'  (92?  -  Wi^  (92?  -  jla22X,  (pi)  ds 

is  a  martingale.  The  continuity  of  this  martingale  derives  from  the  result  given  in  the  following 
lemma  [13]: 

Lemma  4.15  Let  X  be  a  local  martin9ale  such  that  Xo  =  0  and  A  be  a  continuous  increasing 
process  such  that  Aq  =  0.  //exp  {ikJfi  -  (ifc*/2)i4t}  is  a  local  martingale  for  every  k  ^  R+ ,  then  X 
is  continuous  and  At  =  (X}f. 

Thus,  by  (4.37),  we  have  that  W/^  (g2?  —  (g2?  is  a  continuous  submartingale.  This  is  true 

for  every  K  &  R  and,  thus, 

W(92?-'^(9t?  =  (Wt(92)-K?-Wt(g2? 

=  K^-2KWt(g2) 

is  a  continuous  process  for  each  K,  which  implies  that  Wt  (P2)  is  continuous.  ■ 

Lemma  4.14  and  Theorem  1  [9]  immediately  yield 


Corollary  4.16  Wt  is  a  continuous  process  taking  values  in  S'. 


Corollary  4.17  For  T  >  0,  sup«j  ~  ®xp  {-’4"*01  converges  to  zero  in 

probability  for  each  gi  €  S{R^). 

Proof  The  proof  is  completely  analogous  to  that  of  Theorem  9.2.1(b)  [5],  using  the  fact  that,  for 
each  p3  €  S{R'‘),  Wt  (ya)  is  a  continuous  process.  g 


4.5  Uniqueness 

In  this  subsection  we  shall  prove  that  the  martingale  problem  for  A,  introduced  in  Section  3  as: 

(4.38)  ^  =  {  exp{-(ffi,/ii) +  t(ff2,A‘2)}. 

exp  {-  (51, Ml}  +  »■  (52, M2)}  {5  (-Aji  -  2i)igi  +  au5?  -  210,2^152  -  02252)  ,Mi)  : 
ffi€P(iA)+,  «,2€5(ii‘*)}, 

has  a  unique  solution. 

Throughout  this  article  we  worked  with  real-valued  functions.  Now  for  simplicity,  we  switch  to 
complex-valued  functions.  We  obtained  in  Corollary  4.12  that  all  the  limit  points  of  W^”^)} 

are  solutions  of  the  martingale  problem  for  {Ai,  v).  By  the  obvious  fact  that 

{exp{-(pi,/ii)}sin((ff2,M2))U€xp{-(5i,/ii)}cos((p2,M2))  :5i  €  V{A)+,g2  €  5}  C  V{Ai), 

we  obtain  that  all  the  limit  points  are  a.so  solutions  for  (A,  v).  Thus,  in  order  to  prove  weak 
convergence,  it  is  sufficient  to  prove  uniqueness  for  (A,  u).  Let  (X,  W)  be  any  limit  point  of 
{(X^"^,  W^"J)}.  For  simplicity,  in  the  next  two  lemmas  we  shall  use  the  notation  introduced  in 
Section  4.4;  that  is  {g,W{t))  and  {g,X{t))  will  be  denoted  by  Wt{g)  and  Xt{g),  respectively. 

Lemma  4.18  For  each  gi  €  P(5A)+  ,  g2  €  S{R'^),  the  bracket  process  of  Xt  (pi)  and  Wt  {92)  is 

(4.39)  {X{gi),Wig2))t  =  /‘a-. (0125152)  *. 

Jo 


Proof  By  Ito’s  formula 

exp  {-Xt  (si)  -I-  iWt  (S2)}  = 

exp  {-Jfo  (5j )  +  (52)}  +  ®*P  (52)}A^»  (-»7i5i  -  |Affi)  ds 

+  1  /‘exp{-X,(si)  -b  »W.(s2)}d((-A-(5i)  +  iW{g2)),)  -b  M., 
where  Mt  is,  at  least,  a  local  martingale.  Recall  that 
exp  {-Xt  (51)  -b  iWt  (52)} 

-  exp  {-X,  (gi)  +  iW,  {g2))X,  (-2inji  -  A^i  -b  on(ji)’  -  2tai25i52  -  022(52)’)  ds 


is  a  martingale  and 


{Xi9))t  =  J^X,(au9^)  di,  \/g€V{\A), 

(W'Ci;))*  =  £x,{a229^)  ds,  Vff€5(7?‘'). 

Combining  the  above  result  with  the  fact  that 

{-X{gr)  +  tW{g2)),  =  {Xigx)),- (W{g2)),-2i{Xigx),W{g2)),  , 

we  are  done.  ■ 

Let  Co(iZ‘^)  be  the  set  of  continuous  functions  tending  to  zero  at  infinity.  Introduce  the  additional 
notation 

(4.40)  C'{r<^)  =  {9:  9  =  9i  +  i92-9ieC,  (/?-*)  ,  <72  €  Co  (R'^)  }  , 

(4.41)  C'{r<‘)^  =  {5:56  C’{r<^)  ,gieV  (Ia)^  ,  <72  €  D  (^a)  n  Co  (/?“')  }  , 
and  the  operator 

(4.42)  I  =  {  exp {-  (fli,7ii)  +  t  (52, 7^2)}, 

exp {-  +  » (52.7^2)}  (5  (“Aji  -  271^1  +  an9i  “  2ia229i92  -  022p|)  ,7^i)  : 

«7i€C'(ii‘')^,<72€5(ii-')}. 

Then  we  have 

Lemma  4.19  (X,  IV)  is  a  solution  of  martingale  problem  for  (A,  v)  if,  and  only  if,  it  is  a 
solution  of  the  martingale  problem  for  (A,  v). 

Proof  Since  A  C  A,  if  (A^,  IV)  is  a  solution  for  (A,  1/),  then  it  is  also  a  solution  for  (A,  u).  Let 
(A,  IV)  be  a  solution  for  (A,  v).  For  each  91  €  CXR^)+ ,  92  €  we  obtain,  by  Ito’s  formula, 

that 

exp{-A|(ffi)  +  tIVt(52)} 

=  exp  {- Ao  (gi)  +  *IVo  (92)}  +  ^  exp  {-A,  (gi )  +  ilV,  (g2)}X,  (-/h9i  -  j^9i)  ds 
+  i  /‘exp{-A.(si)  +  tIV,(s2)}d((-A(si)  +  *IV(s2)).)  +  Af„ 

=  exp  {-Ao  (ffi)  +  iIVo  (ff2)}  +  5  exp  {-A,  (^i)  +  ilV.  (92)} 

X  X,  (-2t7iJi  -  Ag\  f  011(51)*  —  210125153  -  032(52)’)  ds  +  Mj, 
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where  the  second  equality  follows  by  previous  lemma,  and  Mi  is,  at  least,  a  local  martingale.  But 
all  the  terms  in  the  right  and  left  hand  sides  are  bounded,  so  that  Mt  is  a  martingale.  g 

In  the  proof  of  uniqueness  the  following  lemma  will  play  a  crucial  role: 

Lemma  4.20  The  nonlinear  evolution  equation 


(4.43) 


dU(t) 

dt 

.  m 


-fanf/(t)*  +  £/(t)qi  +  ^At/(t)  +  toi2f/(t)p2  +  §022(52)* 


5i  > 


where  g\  €  CXK‘)+  ,  52  €  5(il^),  has  a  unique  strong  solution  on  R+  ,  such  that  t/(t),  dU  (t)  /dt, 
and  lAU{t)  are  continuous  functions  from  R+  to  C(J2''),  and  U{t)  6  Ct /?'')+  for  each  t  e  R+ . 


Proof  Appendix. 

Lemma  4.21  If  (A’(t),lV(t))  is  a  solution  of  the  martingale  problem  for  (A,  u)  in  i?Afjpx5'[0,c»), 
then,  for  each  gi  €  T  >  0, 

(4.44)  exp  {-  {UiT  -  t),  X(t)>  +  i  (52,  lV(t))} 

is  a  martingale  for  0  <  t  <  T  ,  where  U{t)  is  given  by  the  unique  solution  of  (I. IS)  with 

V{Q)  =  gieC\R<^)^. 


Proof  Define 


(^•^)  ^(91^92)  =  -  §011  (51)*  +  10125152  +  5O22  (52)* . 

Let  A(t)  =  (A(t),W'(t))  be  a  solution  of  the  martingale  problem  for  (A,  1/),  v  €  V{Mf  x  S'). 
Define 

(4.46)  tt  (s,  A(t))  =  exp  {-  {U(T  -  s),X(t))  +  i  {gt,  H'(t))},  VO  <  s  <  t  <  T. 

For  each  {(11,(12)  €  Mp  x  S'  ws  have 


(4.47) 


ds 


-(t;(r-s),pi)  +  i(p2,M2)}^^|_  ^  ,(iiy 


Hence 


^[«(t2,.^(t2))-«(tl,A(t2))|.F„ 


-  Ey^'"eKp{-{U{T-s),X{t2))  +  i{g2,W{t2))}(^^^\^^^  ^  ,^(*2)^  ds| 

=  Eyy{s,X{t2))ds\j^ti'\. 


Fix  0  <  ti  <  r.  Since  gi  €  CXR%,  we  have  U{T-ti)  €  C'(H^)+,  for  every  0  <  ti  <  T.  Therefore, 
by  the  definition  of  the  martingale  problem  for  (A,  i/),  we  have  that 

exp  {-  (U{T  -  h)gi,X{t))  +  i(52,W^(t))} 

-/«*?{-  iU{T  -  ti),Xis))  +  i(52,W'(«)»  X  (-F{U{T  -  *0,52)  -  \AU{T  -  U),X{s))  ds 
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16  a  martingale  for  0  <  t  <  T.  Define  the  function  tr  by  the  following: 

(4.48)  E  [u  (ti,^(t2))  -  ti  (t,.^(«,))|jF,.] 

=  E  (exp  {-  {V{T  -  U),X{t2))  +  i  (p2,  IV(tj))} 

-exp{-  mr  -  ti),X(ti))  +  t{92,W(U))}\J^tr] 

=  E  [j[‘*  exp  {-  {UiT  -  tr),  X^s))  +  i  {g2,W(s))} 

x{-F(UiT-ti),g2)-\AU{T-h),X(s))ds\Et,] 

=  ^  <  U  <  tj  <T  <  oo 

It  is  easy  to  check  that  all  conditions  of  the  Theorem  4.3.4  [5]  are  satisfied  and  hence 

(4.49)  u  (t,X(t))  -  {«  («.^(«))  +  (s,a,^(a))}  ds 

is  an  -martingale.  Finally,  substituting  the  definitions  of  u,  v,  w  we  obtain  the  desired  result.  | 

Theorem  4.22  The  martingale  problem  for  (A,  i/)  has  a  unique  solution. 

Proof  Let  (Jf(t),  W'’(t))  be  any  solution  of  the  martingale  problem  for  {A,v),  v  €  V{Mf  x  S'). 
Take  g^  €  S(R^)  and  51  €  V{^A)+  .  (Recall  this  means  that  gi  =  gn  +  »Pi2  €  C\R'^)+  ,  gu  =  0.) 
Then  by  the  previous  lemma,  setting  T  =  t,  we  obtain 

f^l€xp{-  (51  ,X{t))  +  i{g2,Wm]  =  -Plexp{-  mt), XiO))  -1-  t(52,W^(0))}],  VO  <  t  <  00, 

where  U{t)  is  the  unique  solution  of  (4.43)  with  1/(0)  =  pi.  Thus,  by  Corollary  1.9,  we  obtmn 
that  that  any  two  solutions  of  the  martingale  problem  for  (A,  1/)  have  the  same  one-dimensional 
distributions.  By  Theorem  1.4  the  desired  result  follows.  ■ 

Remark  4.23  The  proof  of  Theorem  S.S  is  now  finished,  since,  by  Corollary  4. 12  and  Theo¬ 
rem  all  the  conditions  of  Theorem  1.3  hold. 


Appendix:  The  solution  of  a  particular  non-linear  evo¬ 
lution  equation 


Lemma  A.l  The  nonlinear  evolution  equation 


(») 


du{t) 


I  m 


— jOu{/(/)*  +  U{t)rfr  +  |A(/(t)  -f  iai2U(t)g2  +  5022(52)* 


ffi. 


where  51  €:CXR'‘)+,  g2  €  S(R^),  has  a  unique  strong  solution  on  R+ ,  such  that  U(t),  dU(t)/dt,  ^AU(t) 
are  continuous  functions  from  R-^  to  C\R^),  and  V{t)  €  CXR^)+  for  each  t  ^  R+  . 
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Recall  that  C\It^)+  is  defined  by  (4.41). 

Lemma  A.2  For  each  gi  €  CXR^)+  there  exists  a  unique  strong  solution  of  (w)  on  (0,  (mar), 
V)here  (m«»  <  »•  Moreover,  if  t^^ax  <  oo,  then  lim  ||l/(()||  =  oo. 


Proof  The  proof  was  outlined  by  Perkins  [12]  for  the  case  of  real-valued  functions.  In  our  case 
the  proof  is  completely  the  same  (c.f.  Theorems  (3.13  —  3.15)  from  [12]).  The  more  complicated 
part  of  existence  of  a  solution  without  explosion  will  be  given  here.  That  is,  we  shall  show  that 
(max  —  OO- 

Before  starting  our  proofs  we  introduce  several  lemmas  from  [12]. 

Lemma  A.S  Let  X”'  be  the  super- Broumian  motion  given  by  the  unique  solution  of  martingale 
problem  for  (A,  dm) ,  m  €  Mp,  where 

A=  |exp{-(5,/i)},  exp  {-  {g,p)}  (-ag  +  \  {txgf  -  \Ag,p.)  :  g  e  V  ^  ,  a,«T€  r|. 

Then  for  each  V> :  [0,oo)  x  — ►  R,  such  that  0(s),  dij}{s)/ds,  ^Ail>{s)  are  strongly  continuotu 

functions  from  [0;  T]  to  Ci{R^)  we  have 

(A.1)  xr  im)  =  mim)  +  +  iA,^.(s)  aV>(s))  ds  ^  Zt  W , 

where  Zt{i{))  is  the  martingale  with  increasing  process 

{z{ii,)),= jy  XT  {ii>{s?)ds. 

Lemma  A.4  The  unique  strong  solution  of  the  nonlinear  evolution  equation 

=  -^<7M0*  +  o^(0  +  |AV(0 

=  9^'Di\A)+ 
satisfies  V(t)  €  Ci{R^)^  for  each  t  >  0. 

Let 


dV(t) 


V(0) 


(A.2) 


m  =  Ui(t)  +  i[/2(i) 


be  the  unique  strong  solution  of  (w)  on  [0,  tmu),  where  Ui  and  I/2  are,  respectively,  the  real  and 
ima^nary  parts  of  the  solution.  We  can  see  from  (w)  that  that  Ui  and  U2  satisfy 


(A.3) 


dt 


-\atiUi(ty  +  Ui{t)in  +  \AUi{t) 


+  Joiiffj(()*  -  ai2V2{t)g2  +  |oj2  (52)* 


ViiQ)  =  gi\  y 


(A.4) 


dU2it) 

at 


—oiiUi{t)U2it)  +  U2{t)rii  -I-  ^AU2{t)  -f  ai2Ui{t)g2 


£^2(0)  =  g\2  • 
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Let  us  define  for  each  h,f€ 

(A.5)  -  ai2h/  +  ^022/^- 


It  is  easy  to  check  that 


FihJ)ix)  >  0,  V*  €  R*. 


Lemma  A.5  Under  the  definitiona  (A.S),  (A.4)  Ui{t)  €  Ci{R^)+,  for  all  t  <  t^ax- 

Proof  We  shall  prove  the  result  by  contradiction.  Assume  that  there  exists  a  t  <  tmaz,  such  that 
Uit,x)  <  0-  Define 


t’ =  inf  |t  <  fmax  :  inf  t^(f,*)<o|. 
I  xeR*  > 


The  solution  {/(•)  of  (x)  is  continuous  on  {0,  tmax),  nnd  pi  €  .  Thus  t“  >  0,  and,  for  each 

s<f. 


(A.8)  f^,(5)€C,(jr')^. 

Define  ^^(s)  =  l/i(f*  —  s).  Obviously  we  have 

^V>(*)  _  _  ^t^i(n)|  w  ,  <  <* 
ds  du  L-f_,  ’  “ 

Now  let  X”*  be  the  super- Brownian  motion  which  is  the  unique  solution  of  the  martingale  problem 
for  (A*,d„),  m  €  Mf,  where 

A’  =  |exp{-  exp{-  (Pii)*  -  :  9u  €  D  (2^)^}  • 

By  Lemma  A.S,  and  equations  (A.3),  (A.5)  we  obtain,  for  all  t  <  t’, 


xrmr-t)) 


=  m(Ui(n)  +  J^XT  -  s)  -I-  rnUi{f  -  s)j  ds  -I-  Z.  ( 

=  m(£^i(r))  +  J*XT  ds  -  jf'jT,”*  (FmT  -  s),^,))  ds  +  Z*  (^) , 


where  Zt  {if>)  is  the  martingale  on  [0,f*]  snth  increasing  process 


{Zii>)),  =  £at2Xr{Ui{f-s)^)  ds. 


Rearanpng,  we  obtain 


-xr  mr  -  0)  -  f  XT  (Fmr  -  S),P2))  ds  =  -m(Ui{f))  -  Zt  (0)  -  \  (Z  (0)). 

JO 


(A.9)  exp  {-Jfr  mr  -  0)  -  j^XT  {Fmt-  -  a),,,)) 

=  exp  {-m{Ui{r)))  exp  {-Zj  (0)  -\{Z  (V’)),}- 

Note  that  exp  {~Zt  (V*)  —  ^{2  (V*))(}  at  least,  a  local  martingale  on  [0,  t*],  but  that  the  left  hand 
side  of  expression  (A.9)  is,  by  (A.6),  (A.8),  bounded  and  so  is,  in  fact,  a  martingale.  Setting  t  =  t* 
we  obtain 

(A.IO)  £;[exp|-^,7(pn)-j[‘^r(TO(t*--),lh))d-|]  =  exp  {-m(f/,(r))}. 

On  the  other  hand,  by  (A.6), 

(A.ll)  E  [exp  j-A,?  (j„)  -  J*  A"  (F(t^,(f  -  s),5,))  dsjj  <  E  [exp  {-A',’?  (p„)  ds}] . 
Define  V{t)  to  be  the  unique  strong  solution  of 

=  -\anV{tf  +  r,iV{t)+\AV{t) 

\  V(0)  =  p„6D(|a)^. 

By  Lemma  A.4  we  have  that  V(t)  €  Ct{R^)+  for  each  t  >  0.  By  standard  arguments  from  the 
theory  of  superprocess  we  find 

(A.12)  £^[exp{-A’,’?(5u)}]  =  exp{-m(V(t*))}. 

Combining  (A.IO),  (A.ll),  (A.12)  we  obtain  that 

exp  {-m(f/i(r))}  <  exp{-m(V(t*))}. 

The  above  expression  holds  for  all  m  €  Mp.  Therefore,  setting  m  =  S^,  we  obtain 


and  hence 


exp{-t^»(<".®)}  <exp{-V^(<*.x)}.  Vx  €  ii*' 


iirf  t^,{t*,x)  >  inf  K(t-,x). 
*€/l^  «€IH 


But  infs£|{a  V{t*,x)  >  0,  so  that  inf^c/ta  (7i(t*,x)  >  0,  which  is  a  contradiction. 
Lemma  A.6  There  exist  constants  Ci ,  Ca  >  0  such  that 


(A.13) 


||t^2(0ll  ^  «P  V t  <  t„ 


and  hence,  if  tmas  <  oo,  then 
(A.14)  1 


limsupllf/aCOii  ^  C’lexplCatm.x}- 


Proof  Consider  two  cases 

(1)  ih  >  0 

Define 

(A.15)  t/2(t)  =  *  <  *moz  • 

By  (A.4)  we  obtain 


dUijt) 

dt 


(—auUi(t)U2it)  +  U2(t)Tfi  +  ^AU2(t)  +  012^^1(0^2)  ®*P  {“’?i0 


-  »?i  1/2(0  exp  {-•nO 


=  —oiiU\(t)U2(t)  +  |A£/2(t)  +  «i2t/i(0ff2exp  {— r/iO 


=  i/l(0  (“®11^2(0  +  0l252exp  +  5At/2(0*  i  <  ^maz  ■ 


Define 

Cl  =  2  max 

Assume  that  there  exists  to  <  tmaz,  such  that 


supi|!/2(o||>C’i. 


Choose  t*  such  that 

(A.16)  ||i/2(C)||  =  sup  l^wl- 

By  choice  of  Ci,  we  have  that  t*  >  0.  For  all  t  <  tmaz>  U2{t)  €  so  that  there  exists  a 

X*  €  R'^,  for  which 

p2(t%**)|=l£/2(r)|. 

Assume,  without  loss  of  generality,  that  x*  is  the  point  of  maximum  of  {/2(t*)  (in  the  case  of  a 
minimum  the  proof  is  analogous).  By  the  previous  lemma,  Ui{t)  €  Ci{R^)+,  for  all  t  <  tmaz-  By 
the  choice  of  Ci  we  obtain  that 

£/i(t*,x*)  (-onC2(f,x*)  +  01292 exp  <  0. 

Since  x*  is  the  point  of  maximum  of  i/2(**)»  the  positive  maximum  principle  ([5],  p.l65)  implies 
that 

iA{/2(f,**)  <  0. 


Hence 


<  0. 


gf/a(f.«*) 

at- 
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Recall  that  1*  >  0,  so  that  there  exists  a  t  <  t*  for  which  U2{i,x’)  >  U2{t',x’).  But  this  contra¬ 
dicts  (A.  16).  Hence  we  obtain  that 

^sup  |t^2(0|  <  Cu 

and  by  the  definition  of  U2  we  have 

(A.17)  ||t^2(0ll  <  Cl  exp  t  <  t„ax , 

which  was  what  we  wanted  to  prove. 

(2)  fh  <  0 

The  proof  is  the  same,  the  only  difference  being  that  we  need  not  introduce  (/j,  as  in  the  previous 
case.  ■ 


Proof  of  Lemma  A.l  Assume  by  contradiction  that  (max  <  00.  It  is  well-known  [11]  that  if 
(/(t)  is  a  strong  solution  of  (w)  on  [0,  tmax)  then  it  is  also  a  mild  solution  of  the  following  integral 
equation  : 

[/(i)  =  5(l)fli  +  ~  +  C(s)t}i  -f  iai2U(s)g2  +  ^022(93)^)  da, 

where  S(t)  is  the  semigroup  generated  by  |A.  By  (A.2),  (A.3)  we  obtain 
(A.18)  Uiit)  =  Sit)gu  J^Sit  -  a)  {-ia„t;,(s)*  +  Ui{a)rjt  +  F(f/2(s),pj))  da 
Recall  that,  by  Lemma  A.5,  we  have  Ui{t)  €  C\R'^)+  for  all  t  <  tmax-  Thus 
0  <  Ui{t)  =  5(t)pn  -f  S{t  —  s)  da  +  S{t  -  a)F  (1^2(5)15^2)  da 

<  Sit)gu  +  f  Sit-  s)iUtia)ni)  da+  f  S{t  -  s)F  (1^2(5),  52)  da 

Jo  Jo 

<  llffiill  +  Iml  /  IICi(s)ll  da -hi  sup  IIF(1/2(s),  92)11 

Jo  »<f 

Finally,  defining  C  =  fma»8up,<,^  Il^(f^2(5)>ff2)||.  (C  <  00  by  Lemma  A.6)  and  A  =  ||pii||  +  C 
we  obtain 

(A.19)  \\Uiit)\\  <  K  +  \fn\  /‘  ||t^i(5)||  da,  Vt  <  Uax , 

Jo 

mnd  by  Gronwall^s  inequality  we  have  that 

||I^l(01l  <  Aexp  {|;;i  1 1},  Vt  <  tmax. 


This  implies  that 

(A.20)  limsup||£fi(t)||  <  AexpdiTiltmxx}  <  00. 
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Combining  the  above  result  with  Lemma  A.6  we  obtain 


* 


(A.21)  lira  ||f/^(t)ll  5  limsup ||l/i(t)||  +  limsup ||l72(t)||  <  00, 

‘T*""  «ttmM  ItfmM 

which  contradicts  Lemma  A.2.  Hence  tmax  =  <X),  which  finishes  the  proof  of  the  lemma. 

■ 
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